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Abstract 

This overview is concerned with the well-posedness problem for the isentropic com¬ 
pressible Euler equations of gas dynamics. The results we present are in line with the pro¬ 
gram of investigating the efficiency of different selection criteria proposed in the literature 
in order to weed out non-physical solutions to more-dimensional systems of conservation 
laws and they build upon the method of convex integration developed by De Lellis and 
Szekelyhidi for the incompressible Euler equations. Mainly following [5], we investigate 
the role of the maximal dissipation criterion proposed by Dafermos in [6]: we prove how, 
for specific pressure laws, some non-standard (i.e. constructed via convex integration 
methods) solutions to the Riemann problem for the isentropic Euler system in two space 
dimensions have greater energy dissipation rate than the classical self-similar solution em¬ 
anating from the same Riemann data. We therefore show that the maximal dissipation 
criterion proposed by Dafermos does not favour in general the self-similar solutions. 


1 Introduction 

We consider the isentropic compressible Euler system of gas dynamics in 2 space dimensions 
(cf. [8] or [15] or [1]). It is obtained as a simplification of the full compressible Euler equations, 
by assuming the entropy to be constant. The state of the gas is described through the state 
vector 

V = ( p,v ) 

whose components are the density p and the velocity v. The system consists of 3 equations 
which correspond to balance statements for mass and linear momentum. The corresponding 
Cauchy problems reads as 

p + div x (pv ) = 0 

(pv) + div^ (pv <g> v) + V x [p(p)] = 0 

•,0) = p° 

. .-,0) = v°, 
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with t G R + , x G M 2 . The pressure p is a function of the density p determined from the 
constitutive thermodynamic relations of the gas under consideration and it is assumed to 
satisfy p' > 0 (under this assumption the system is hyperbolic). We will work with pressure 
laws p(p) = p 1 with constant 7 > 1 . 

Our aim is to discuss the issue of uniqueness of weak solutions to the Cauchy problem 
(1.1). The theory of the Cauchy problem for hyperbolic systems of conservation laws is 
typically confronted with two major challenges. First, it is well-known that classical solutions 
develop discontinuities, even starting out from smooth initial data. In the literature this 
behaviour is known as breakdown of classical solutions. Therefore, it becomes imperative to 
introduce the notion of weak solution. However, weak solutions fail to be unique. In order to 
restore uniqueness restrictions need to be imposed in hope of singling out a unique physical 
solution. When dealing with systems coming from Physics, as in our case, the second law of 
Thermodynamics naturally induces such restrictions, such admissibility criteria by stipulating 
that weak solutions are admissible/entropy solutions if they satisfy some entropy inequalities 
(see (2.3) for the specific case of the compressible Euler system). Finally, a third important 
challenge then arises: do entropy inequalities really serve as selection criteria? Are admissible 
solutions unique? Or at least, do there exist efficient criteria restoring uniqueness? This is 
a central problem to set down a complete theory for the Cauchy problem. It has deserved 
a lot of attention, but positive answers were found only for scalar conservation laws or for 
systems in one space dimensions (under smallness assumptions on the initial data): for a 
complete account of the existing literature we refer the reader to [ 8 ] and [15]. When dealing 
with systems of conservation laws in more than one space dimension, it is still an intriguing 
mathematical problem to develop a theory of well-posedness for the Cauchy problem which 
includes the formation and evolution of shock waves. 

In 2006, Elling [12] studied numerically a particular case of initial data for the two dimen¬ 
sional non-isentropic Euler equations. His results show that the numerical method does not 
always converge to the physical solution. Moreover, they suggest that entropy solutions (in 
the weak entropy inequality sense) to the multi-dimensional Euler equations are not always 
unique. 

In a groundbreaking paper [10], De Lellis-Szekelyhidi give an example in favour of the 
conjecture that entropy/admissible solutions to the multi-dimensional compressible Euler 
equations are in general not unique. The non-uniqueness result by De Lellis-Szekelyhidi is 
a byproduct of their analysis of the incompressible Euler equations based on its formulation 
as a differential inclusion (see [9] and [11]) combined with convex integration methods: they 
exploit the result for the incompressible Euler equations to exhibit bounded initial density and 
bounded compactly supported initial velocity for which admissible solutions of ( 1 . 1 ) are not 
unique (in more than one space dimension). However the initial data constructed in [10] are 
very irregular. The result by De Lellis and Szekelyhidi is improved by the author in [2] where 
it is proven that non-uniqueness still holds in the case of regular initial density (see also [4] for 
further generalizations). Non-unique solutions constructed via convex integration are referred 
to as non-standard or oscillatory solutions. Moreover in [3], using the Riemann problem as 
a building block, the authors show that, in the two dimensional case, the entropy inequality 
(see (2.3)) does not single out unique weak solutions even under very strong assumptions on 
the initial data ((p°,u°) G VE 1 , 00 (M 2 )): 

Theorem 1.1 (Chiodaroli, De Lellis, Kreml) There are Lipschitz initial data (p°,v°) 
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for which there are infinitely many bounded admissible solutions (p,v) of ( 1 . 1 ) on M 2 x [0, oo[ 
with inf p > 0. These solutions are all locally Lipschitz on a finite interval on which they all 
coincide with the unique classical solution. 

This is proven by constructing infinitely many entropy weak solutions in forward time to a 
Riemann problem for (1.1) whose Riemann data can be generated, backwards in time, by a 
classical compression wave: the Lipschitz initial data of Theorem 1.1 will be provided by the 
values of the compression wave at some finite negative time. It is clear now that the infinitely 
many admissible solutions constructed in Theorem 1.1 all coincide with the unique classical 
solution (compression wave) on a finite time interval whereas non-uniqueness arises after the 
first blow-up time. 

This series of negative results concerning the entropy inequality as selection criterion for 
system ( 1 . 1 ) motivated the authors to explore other admissibility criteria which could work 
in favour of uniqueness, in particular we investigated an alternative criterion which has been 
proposed by Dafermos in [ 6 ] under the name entropy rate admissibility criterion. The ideas 
developed in [3] enabled us to prove in [5] the following theorem: 

Theorem 1.2 (Chiodaroli, Kreml) Let p(p) = p 1 , 1 < 7 < 3. There exist Riemann data 
for which the self-similar solution to ( 1 . 1 ) emanating from these data is not entropy rate 
admissible. 

This result does not exclude that the entropy rate admissibility criterion could still select a 
unique solution, but surely prevents the self-similar solution to be the selected one. Moreover, 
since Theorem 1.2 is proven using non-standard solutions as competitors, with respect to 
Dafermos’ criterion, for the self-similar solution, we can affirm that the entropy rate criterion 
cannot, at least in our setting, rule out oscillatory solutions obtained via convex integration. 
In the rest of the paper we will further explain this result. 


2 Entropy criteria 

2.1 Entropy inequality 

We recall here the usual definitions of weak and admissible solutions to (1.1) in the two- 
dimensional case. 


Definition 2.1 By a weak solution of (1.1) on M 2 x [0, 00 [ we mean a pair ( p , v) £ L°°(R 2 x 
[0, oo[) such that the following identities hold for every test functions £ (^(M 2 x [0, oo[), 
£ C c °°(M 2 x [0,oo[); 


( 2 . 1 ) 


/ 0 Jr 2 


[pdt'fi + pv • V x ip\ dxdt + / p°(x)'fi(x,0)dx = 0 

Jr 2 


poo P P 

(2.2) / / \pv ■ dtfi + pv ® v : D x (f> + p(p)div x (f>\ + / p°(x)v°(x) ■ fi(x,0)dx = 0. 

J 0 Jr 2 Jr 2 
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Definition 2.2 A bounded weak solution (p,v) of (1.1) is admissible if it satisfies the fol¬ 
lowing inequality for every nonnegative test function tp G C£° (M 2 x [0, oo[): 


/ o Jr 2 


Mp) + P~w~ \ d hP + pe(p) + p-rr +p{p) v • V x p 


(2.3) 


+ 


p 0 {x)e(p°(x)) + 


) ^r 


<p(x, 0) dx > 0 . 


Note that (2.3) is rather a weak form of energy balance. 


2.2 Entropy rate admissibility criterion 

An alternative criterion to the entropy inequality has been proposed by Dafermos in [6] under 
the name of entropy rate admissibility criterion. In order to formulate this criterion for the 
specific system (1.1) we define the total energy of the solutions (p, v) to (1.1) as 

(2.4) E\p, u](t) = J ^ (pe{p) + dx. 


Let us remark that in Dafermos’ terminology E[p,v\(t) is called “total entropy” (see [6]). 
However, since in the context of system (1.1) the physical energy plays the role of the math¬ 
ematical entropy, it is more convenient to call E[p,v](t) total energy. The right derivative of 
E[p,v\(t) defines the energy dissipation rate of (p,v) at time t: 

(2.5) D[M(<) ^ d+E j’” ■ 

Since our solutions will have piecewise constant values of p and |u| and it is easy to see that 
the total energy of any solution we construct is infinite, we shall restrict the infinite domain 
M 2 to a finite box (— L, L) 2 and denote 

(2.6) E L \p,v]{t) = ^ ^ bsfp) + pL-J 

(2.7) o 1 [p,„]( t) = d ±^LrM. 

The problem of infinite energy of solutions may be solved also by restricting to a periodic 
domain and constructing (locally in time) periodic solutions. This procedure is carefully 
described in [5]. 

Using the concept of energy dissipation rate, Dafermos in [6] introduces a new selection 
criterion for weak solutions which goes under the name of entropy rate admissibility criterion. 
We recall here the definition of entropy rate admissible solutions. 



Definition 2.3 (Entropy rate admissible solution) A weak solution (p,v) of (1.1) is 
entropy rate admissible if there exists L* > 0 such that there is no other weak solution 
(p,v) with the property that for some r > 0, (p,v)(x,t) = (p,v)(x,t) on M 2 x [0,r] and 
D l [p,v\(t ) < D l [p,v\(t) for all L > L*. 

In other words, we call entropy rate admissible the solution(s) dissipating most total energy. 
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3 Background literature and main results 


The investigation of the entropy rate admissibility criterion initiated with the paper [ 6 ] of 
Dafermos where he puts it forward and moreover proves that for a single equation the entropy 
rate criterion is equivalent to the viscosity criterion in the class of piecewise smooth solutions. 
Later on, following the approach of Dafermos, Hsiao in [14] proves, in the class of piecewise 
smooth solutions, the equivalence of the entropy rate criterion and the viscosity criterion 
for the one-dimensional system of equations of nonisentropic gas dynamics in lagrangian 
formulation with pressure laws p(p) = p 7 for 7 > 5/3 while the same equivalence is disproved 
for 7 < 5/3. For further analysis on the relation between entropy rate minimization and 
admissibility of solutions for a more general class of evolutionary equations we refer to [7]. 
However, to our knowledge, up to some time ago the entropy rate criterion had not been tested 
in the case of several space variables and on broader class of solutions than the piecewise 
smooth ones. 

Recently Feireisl in [13] extended the result of Chiodaroli [2] and obtained infinitely many 
global admissible weak solutions of ( 1 . 1 ) none of which is entropy rate admissible: this results 
seems in favour of the effectiveness of the entropy rate criterion to rule out non-standard 
solutions (i.e. constructed by the method of De Lellis and Szekelyhidi). In [5] we have 
actually shown that for specific initial data, and in the two-dimensional case, the oscillatory 
(non-standard) solutions dissipate more energy than the self-similar solution which may be 
believed to be the physical one. The results obtained in [3] hinge upon some of the ideas 
devised in [3] combined with novel developments to deal with the entropy rate criterion. 

We refer also to the work [16], where Szekelyhidi constructed irregular solutions of the 
incompressible Euler equations with vortex-sheet initial data and computed their dissipation 
rate. 

We focus on the Riemann problem for the system (1.1),(2.3) in two-space dimensions. 
Hence, we denote the space variable as x = (aq, x 2 ) £ M 2 and consider initial data in the form 

( (, p-,v _) if x 2 < 0 

(3.1) (j>°(x),v°(x)) := 

{ (p+,v+) if x 2 > 0 , 

where p±,v±i,v± 2 are constants. Our concern has been to compare the energy dissipation 
rate of standard self-similar solutions associated to the Riemann problem (1.1), (2.3), (3.1) 
with the energy dissipation rate of non-standard solutions for the same problem obtained by 
the method developed in [3]. 

We obtained the following results. 

Theorem 3.1 Let p(p) = p 1 with 7 > 1. For every Riemann data (3.1) such that the self¬ 
similar solution to the Riemann problem (l.l)-(2.3), (3.1) consists of an admissible l—shock 
and an admissible 3—shock, i.e. V-i = v + \ and 



there exist infinitely many admissible solutions to (l.l)-(2.3), (3.1). 
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Theorem 3.1 can be viewed as an extension of the results obtained together with De Lellis 
in [3]. As a consequence of Theorem 3.1 and by a suitable choice of initial data, we can prove 
the following main theorem. 

Theorem 3.2 Let p(p) = p 1 , 1 < 7 < 3. There exist Riemann data (3.1) for which the 
self-similar solution to ( 1 . 1 ),( 2 .3) emanating from these data is not entropy rate admissible. 

Theorem 3.2 ensures that for 1 < 7 < 3 there exist initial Riemann data (3.1) for which 
some of the infinitely many nonstandard solutions constructed as in Theorem 3.1 dissipate 
more energy than the self-similar solution, suggesting in particular that the Dafermos entropy 
rate admissibility criterion would not pick the self-similar solution as the admissible one. 

4 Strategies of proof 

In this section we explain how to prove Theorem 3.1 and 3.2. For the complete proofs we 
refer the reader to [5]. 

Both theorems stem from the framework introduced in [3] where the Riemann problem 
constitutes the starting point for constructing non-unique admissible non-standard solutions. 
In particular, in [3], the authors jointly with Camillo De Lellis developed a method which 
allows to obtain infinitely many entropy (oscillatory) solutions to a Riemann problem provided 
a suitable admissible subsolution exists. 

4.1 From subsolutions to solutions 

We shall introduce the notion of fan subsolution and admissible fan subsolution as in [3, 
Section 3]. 

Definition 4.1 (Fan partition) A fan partition 0 /M 2 x (0, 00 ) consists of three open sets 
P-,Pi, T+ of the following form 


(4.1) 


= (OM) 

t > 0 

and 

X2 < V-t} 

(4.2) 

Pi 

= {(*,*) 

t > 0 

and 

V-t < X2 < v+t} 

(4.3) 

P+ 

= {(*,*) 

0 

A 

•HO 

and 

X2 > n + t}, 


where i/_ < is an arbitrary couple of real numbers. 

We denote by c?Q y2 the set of all symmetric 2x2 matrices with zero trace. 

Definition 4.2 (Fan subsolution) A fan subsolution to the compressible Euler equations 
(1.1) with initial data (3.1) is a triple ( p,v,u) : M 2 x (0,oo) —> (M + ,M 2 ,5q X2 ) of piecewise 
constant functions satisfying the following requirements. 

(i) There is a fan partition P-,P\,P + 0 /M 2 x (0,oo) such that 

(p,v,u) = (p-,V-,U-)lp_ + (pi,V 1 ,Ul)l Pl + (p + ,v + ,u+)l P+ 

where pi,v\,u\ are constants with pi > 0 and u± = v± < 8 > v± — 51 "U± 1 2 Id; 
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(ii) There exists a positive constant C such that 


(4.4) v\ <S> v\ — u\ < —Id; 

(in) The triple ( p,v,u ) solves the following system in the sense of distributions: 

(4.5) dtp + div x (pv) = 0 

(4.6) d t (pv) + div* (pu) + V x (p{p) + ^ (Cp 1 l Pl + p\v\ 2 l P+uP _)) = 0. 


Definition 4.3 (Admissible fan subsolution) A fan subsolution (p, v, u ) is said to be ad¬ 
missible if it satisfies the following inequality in the sense of distributions 


d t {pe{p)) + div x \{pe{p) +p(p))v\ + d t < p— 1p +u p_ +div x p—vl PuP _ 


\v \ 2 


\v \ 2 


(4.7) 


+ 


C 


d t 1 pi — lp, + div x pi v — 1 P] 


C 


< 0 . 


The strategy which lies behind Theorem 3.1, as well as behind Theorem 1.1 in [3], consists 
in reducing the proof of the existence of infinitely many admissible solutions to the Riemann 
problem for ( 1 . 1 ) to the proof of the existence of an admissible fan subsolution as defined 
in Definition 4.3. This is the content of the following Proposition which represents the key 
ingredient of [5] and is proven in [3]. 


Proposition 4.1 Let p be any C 1 function and (p±,v±) be such that there exists at least one 
admissible fan subsolution ( p,v,u ) of (1.1) with initial data (3.1). Then there are infinitely 
many bounded admissible solutions {p,v) to (l.l)-(2.3), (3.1) such that p = p and\v\ 2 1 p 1 =C. 

Roughly speaking, the infinitely many bounded admissible solutions (p, v ) of Proposition 
4.1 are constructed by adding to the subsolution solutions to the linearized pressureless incom¬ 
pressible Euler equations supported in Pi. Such solutions are given by the following Lemma, 
cf. [3, Lemma 3.7]. 

Lemma 4.2 Let (v,u) G M 2 x 5q X2 and Co > 0 be such that v <8> v — u < ^Id. For any 
open set Q C R 2 x M there are infinitely many maps ( v,u ) G L°°(M 2 xR,l 2 x Sq X 2 ) with the 
following property 

(i) v and u vanish identically outside Ll; 

(ii) div^u = 0 and dtv + div^u = 0; 

(in) (v + v) (g> (v + v) — (u + u) = ^Id a.e. on D. 

Proposition 4.1 is then proved by applying Lemma 4.2 with Ll = P\, ( v,u ) = (ui,«i) and 
Cq = C. It is then a matter of easy computation to check that each couple (p, v + v) is indeed 
an admissible weak solution to (1.1)—(2.3) with initial data (3.1), for details see [3, Section 
3.3]. The whole proof of Lemma 4.2 can be found in [3, Section 4], 
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4.2 Concluding arguments 

Thanks to Proposition 4.1, Theorem 3.1 amounts to showing the existence of a fan admissible 
subsolution with appropriate initial data under the hypothesis that (3.1) is such that the self¬ 
similar solution to the Riemann problem (1.1), (2.3), (3.1) consists of an admissible 1—shock 
and an admissible 3—shock 

Indeed, a fan admissible sunsolution with initial data (3.1) is defined through a the set of 
identities and inequalities which we recall here (see also [3, Section 5]). 

We introduce the real numbers a, (3, 71,72, u_i, v- 2 , u+i, v +2 such that 

vi = (a,P), 

V- = (v-i,V- 2 ) 

V+ = (v + l,V +2 ) 

ttl = ( 71 72 V 

V 72 -71 / 

Then, Proposition 4.1 translates into the following set of algebraic identities and inequal¬ 
ities. 


(4.8) 

(4.9) 

(4.10) 

(4.11) 


Proposition 4.3 Let P_, Pi, P + be a fan partition as in Definition f.l. 

The constants vi, V-, v+, u\, p_, p + , p\ as in (4.8)-(4.11) define an admissible fan subsolution 
as in Definitions f.2-f.3 if and only if the following identities and inequalities hold: 

• Rankine-Hugoniot conditions on the left interface: 

(4.12) v-(p--pi) = p_v- 2 -pil3 

(4.13) v-(p_v-i - pi a) = p-V-iV-2 - P 172 

c 

(4.14) u_(p_v- 2 - piP) = p-v 2 _ 2 + pi'yi +p{p-) -p(pi) ~ Pi— ; 

• Rankine-Hugoniot conditions on the right interface: 

(4.15) u + {p 1 - p + ) = pip - p + v +2 

(4.16) v+{p\a - p+v+ 1) = P172 - p + v +1 v +2 

c 

(4.17) v+(piP - P+v+ 2 ) = -pi'yi-p+v 2 +2 +p(pi)-p(p+) +pi—\ 


Subsolution condition: 


(4.18) 

(4.19) 


a 2 + (3 2 < C 


C 


C 


— - a + 7! — - P - 71 - (72 - a-fi) > 0 ; 


Admissibility condition on the left interface: 


iy-(p-e(p-) - pie(pi)) + v- p_ 


V- 


C 


P 1- 


1 11 ^ c 

(4.20) < [(p-e(p-) + p(p-))v - 2 - (pie(pi) +p(pi))/3\ + P-V- 2 - pif3— 



Admissibility condition on the right interface: 


u + ( Pl £( Pl ) - p + s{p + )) + v+ 



hf 

2 


(4.21) 


< [(Pl£(pl) +P(Pl))P - {P+Z(P+) + p(p+))v +2 ] + 


(piPfj - P + v +2 



Theorem 3.2 is then a corollary of the following theorem, proven in [5] and showing the 
existence of an admissible fan subsolution, combined with Proposition 4.1. 


Theorem 4.1 Let p(p) = p 7 with 7 > 1. For every Riemann data (3.1) such that V-i = u+i 
and 


(4.22) 


V +2 ~ V-2 < 



P+){P(P-) ~P(P+)) 
P-P+ 


there exist v + , v\, u±, pi, C such that (4.12)-(4.21) hold. 


It remains to prove Theorem 3.2. This is obtained by showing that among the infinitely 
many admissible solutions provided by Theorem 3.1 one has lower energy dissipation rate 
than the self-similar solution emanating from the same Riemann data, thus contradicting 
Definition 2.3. Let us remark that the Riemann data allowing for the result of Theorem 3.2 
are of the same type as the ones of Theorem 3.1, i.e. they admit a forward in time self-similar 
solution consisting of two shocks. We also underline that the self-similar solution depends in 
fact only on one variable, specifically on X 2 - 

Assume from now on for simplicity that V-± = u+i = 0 in (3.1). Let us denote the 
self-similar solution emanating from the Riemann data as in Theorem 3.1. The value of 
the dissipation rate Di[ps,vs\(t) has a specific form for the solution ( ps,vs ) consisting (by 
assumption) of two shocks of speeds 17 and 1/2 • Denoting the middle state (p m ,v m = (0,«)) 
and introducing the notation 

v 2 

(4.23) E± := p±e{p±) + p± y 

(4.24) E m := Pm^iPm) T Pm ~ 
we have 


(4.25) D L [p s ,v s \ = —2L (i 7(.E_ - E m ) + v 2 {E m - E + )). 

Now let us consider a solution (p,v) with the same initial data (3.1) but constructed by the 
method of convex integration starting from an admissible fan subsolution using Proposition 
4.1. We also assume, that the fan admissible subsolution (which exists by Theorem 4.1) has 
an underlying fan partition defined by the speeds z^_ and Although v is not constant 
in Pi we still have, by construction (see Proposition 4.1), that \v\ 2 lp 1 = C, in particular 
Ei := pie(pi) + pi $( is constant in P±. The dissipation rate for all solutions constructed from 
a given subsolution hence depends only on the underlying subsolution and is given by 

(4.26) D L [p,v\ = -2P(z/_(P_ - Pi) + u + (Ei - E + )). 
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If, for a moment, we assume that the speeds of the self-similar solution and of the subsolution 
coincide, i.e. u_ = v\ and v + = z^, it is clearly enough to achieve E\ > E m in order to prove 
Theorem 3.2. Of course, as one can see from [5, Section 4] this is not the case; nonetheless 
the proof works along the same line: we can prove that there is still some freedom in the 
choice of the parameters defining the underlying subsolution for (p, v ) which allows to tune 
them in such a way that Theorem 3.2 holds. For a complete and rigorous proof we refer to 
[5, Section 5]. 
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